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Abstract 

Gamma ray bursts are excellent candidates to constrain physical 
models which break the Lorentz symmetry. We consider deformed dis- 
persion relations which break boost invariance and lead to an energy- 
dependent speed of light. In these models, simultaneously emitted 
photons from cosmological sources reach Earth with a spectral time 
delay that depends on the symmetry breaking scale. We estimate 
the possible bounds which can be obtained by comparing the spec- 
tral time delays with the time resolution of available telescopes. We 
discuss the best strategy to reach the strongest bounds. We compute 
the probability of detecting bursts that improve the current bounds. 
The results are encouraging, depending on the model, it is possible 
to build a detector that within several years will improve the present 
limits of 0.015 m p i. 

1 Introduction 

One of the open questions of high energy physics is how to unify gravity 
with quantum physics. A lot of effort has been devoted to develop a theory 
of quantum gravity. This theory is likely to require a drastic modification 
of our current understanding of the space-time. At present there are two 
formal mathematical approaches : loop quantum gravity and superstring 
theory. Whatever might be the right description of the space-time at very 
short scales, there are some likely physical manifestations. It has been sug- 
gested, for instance, that such a theory would break what we believe to be 
basic symmetries of nature. In [HE], it was shown that Einstein Lagrangian 
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allows for large fluctuations of the metric and the topology of the space- 
time on scales of order of the Planck length, creating a foam-like structure 
at these scales. It has been proposed that the propagation of particles in 
a foamy space-time is strongly affected on short scales. The medium re- 
sponds differently depending on the energy of the particle, in analogy with 
the propagation through a conventional electromagnetic plasma [Hill]. Thus 
space-time might exhibit a non-trivial dispersion relation in vacuum, violat- 
ing therefore Lorentz invariance. 

There are many different ways of breaking Lorentz symmetry; a back- 
ground tensor field, like a magnetic field, breaks the vacuum rotational in- 
variance, for instance. However, it has been shown that there are 46 different 
ways in which the standard model Lagrangian can be modified while remain- 
ing renormalizable, invariant under SU(3) x SU(2) x U(l) and rotationally 
and translationally invariant in a preferred frame 0. Among other effects, 
these terms cause the velocity of light to differ from the maximum attainable 
velocity of a particle, therefore changing the kinematics of particle decays. 
Modifying the dispersion relations of photons and electrons allows for new 
QED vertex interactions like photon splitting in vacuum, vacuum Cerenkov 
effect for electrons, photon decay, electron-positron annihilation to a single 
photon, etc. Ej. 

In this work we consider only rotationally invariant deformations of the 
photon dispersion relation which produce an energy- dependent speed of light. 
If this effect is present in nature, it has to be absent below some energy 
scale, £m p i, (where m v \ is the Planck energy and £ a dimensionless constant), 
high enough to have been undetected so far. Traditionally this energy is 
believed to be the Planck energy which, at first, seems to make hopeless 
any experimental attempt to test these models. However, in 1997 Amelino- 
Camelia et al. [S] suggested that such models can be explored by studying 
the propagation of photons emitted from a distance source like a gamma ray 
burst (GRB) GRBs are short and intense pulses of soft gamma rays that 
arrive from cosmological distances from random directions in the sky. The 
bursts last from a fraction of a second to several hundred seconds. Most 
GRBs are narrowly beamed with typical energies around 10 51 ergs, making 
them comparable to supernovae (for a recent review see [HI CD]). Because of 
the large distances traveled by the photons, these bursts are valuable tools 
to explore energies far beyond the reach of any laboratory on Earth. 

If one considers photons with energies much smaller than the symmetry 
breaking scale, it is possible to expand the dispersion relation in powers of 
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E/^nipi. The first correction produces a tiny departure from the Lorentz 
invariant (LI) equations and we expect that the low energy limit of the de- 
formed dispersion relation can be written generically as: 



E 



) 



n 



E 2 -p 2 c 2 ~eE 2 



(1) 



where e = ±1 takes into account the possibility of having either infraluminal 
or superluminal motion, the latter appearing in some models of quantum loop 
gravity [TTl I12j . Photons simultaneously emitted from a GRB with different 
energies will travel at different speeds, and therefore will show on Earth a 
time delay. The goal of this paper is to explore these high energy Lorentz 
violation models by studying such time delays. We analyze the potential of 
detecting observational consequences of a modified dispersion relation like 



The paper is organized as follows; in section[2]we review cosmological pho- 
ton propagation in the LI theory and show how these results are modified 
when a non LI term is introduced. We compute the travel time of cosmo- 
logical photons and, as a check, compare it with the travel time obtained in 
the Newtonian approximation. In section |3] we turn to the observations and 
show how such models can be tested using GRB observations. We compare 
our method with previous works in section |S] and finally conclude in section 



2 Propagation of the photons 

We consider first the propagation of a particle in a FRW universe, described 
by the metric ds 2 = —c 2 dt 2 + a(t) 2 dx 2 . The Hamiltonian of a relativistic 
particle is 



where p is the comoving momentum and m the mass. The Hamiltonian 
depends explicitly on t through a(t), expressing the fact that the momentum 
is redshifted due to the cosmological expansion. The trajectory of the particle 
is 



Eq. □ 




(2) 




p = constant . 



(3) 
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For a massless particle, like a photon, Eq. El becomes 



ax = c . (4) 

Hence the speed of photons is an universal constant, c, which does not depend 
on the energy. 

When Lorentz symmetry is broken this result is modified. As long as a 
theory of quantum gravity is not available, the high energy corrections to 
the Hamiltonian defined in Eq. El cannot be calculated. Here we will adopt a 
phenomenological approach, assuming that the Hamiltonian at high energy is 
an unknown function of the momentum, which reduces at low energies to Eq. 
El At small energies compared to the symmetry breaking scale, E -C £m p i, 
a series expansion is applicable. We will consider the first order correction 
to the LI theory. 

We consider models which break boost invariance but keep rotational and 
translational invariance. Inspired by Eq. [2] we therefore postulate 



n 



2 2/ / \ n\ -i 1/2 

2 4 , P C f PC 

m c H — 1 + 



n=l,2,... . (5) 



Note that given Eq. ^ there is some arbitrariness in the choice of Eq. E] 
concerning the a~ n factor. We believe that this choice is physically the best 
motivated because any p dependence should be redshifted due to the cosmic 
expansion. 

We define the parameter \i as the ratio between the energy of the photon 
and the Planck energy, 

fi(a,p) = . (6) 

m p ia 

From Eq. El we deduce the new trajectory of a particle 



x(t,p) 



2 

pc 



l + kWl + n + 



2 r ' * i- - m 2 c 4 + £ !|! 



dt 



p = constant . (7) 

We made a linear expansion with respect to fi n since Eq. is only valid to 
linear order. The linear approximation used in equation El remains valid for 
/i < 1. Notice that in the limit when the symmetry breaking scale goes to 
infinity, £ — > oo, we recover Eq. Eland the Lorentz symmetry is restored. 
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The main and striking difference with respect to the LI theory is that the 
speed of a massless particle depends on its momentum 



ax = c(l + ^(l+n)fi(a,p) n C n ) ■ (8) 
2.1 Time delay 

Because of the energy dependence of the speed of light in Eq. |HJ two photons 
emitted at the same time from the same source with momenta p\ and P2, 
will reach Earth at times t\ and t 2 . The comoving distance traveled by both 
photons is the same 

x{ti,pi) = x(t 2 ,p 2 ) , (9) 
t 2 = h + At. 

Rewriting equation [7| in terms of the redshift and particularizing for photons, 
m = 0, one obtains 

Ho Jo \ 2 J .Q m i + Z y + Q A 



where the cosmo logical parameters Hq, Q m and Q\ are evaluated today ^B] 
and we set ao = 1. Notice that 

fi(a,p)=fio(p)(l + z) , fio(p) = /i(a ,p) . (11) 

When £ — > oo, Eq. ITUl becomes the standard definition of the cosmo logical 
distance, as it should. 

Expanding Eq. |§] for small At, we obtain the time delay between two 
photons with momenta p\ and p 2 , 



AtM ^ i±2. Afi" ^j±£3 r , < 1+ ^ , (i2) 



where A// 1 = [^{p-z) — p^ipi)- Fig. [TJdepicts the time delay as a function of 
the redshift and the momentum of the photon. 
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Figure 1: Curves of constant time delay as a tunction of the redshift of the 
source and the energy of the photon. In all the curves n — 1 and £ = 1. 



2.2 Newtonian approximation 

It is instructive to compare the former results with those obtained in the 
Newtonian approximation. For small redshifts one can neglect the expansion 
of the universe and suppose that energy of the photons is constant. The 
delay between a low energy photon traveling at the standard speed of light 
c and the a high energy photon traveling at the modified speed of light v is 

d d . . 

At , 13 

c v 

where v = dE/dp. In this approximation a linear relation between distance 
and redshift holds 

Hence the time delay is 



n 



At — ■ — /igr w - (is) 

2H \fVt m + f^A 

Comparing Eq. E3and Eq. ^]we verify that the Newtonian analysis agrees 
to first order with its relativistic counterpart, Eq. 
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3 Observational Detection of Lorentz Viola- 
tion 



To obtain an experimental bound on the symmetry breaking scale £, we need 
to compare the delay produced by the modified speed of light, Eq. |S1 with 
the time resolution of the observing telescope. For a successful detection of a 
Lorentz violation, the delay has to be larger than the time resolution of the 
telescope: 

At dcl > At ms . (16) 



3.1 Telescope time resolution 

The time resolution of a telescope depends on two factors. The first is the 
intrinsic detector minimal time resolution, Atdetector; which is typically of 
order 10~ 3 — 10~ 4 s. The second factor, At res , is inversely proportional to the 
photon detection rate which, in turn, depends on the detector effective area 
and on the luminosity, spectrum and distance of the source, 

M - = A^kM ' (17) 

where P 7 (Ei, E 2 ) is the photon peak flux in the energy band (Ei, E 2 ) and A 
is the detector effective area. The factor b takes into account the minimum 
number of photons needed to resolve a peak. If the noise of the detector 
is negligible, b is of order 5-10. In the following examples we consider an 
idealized detector with no noise and set 6 = 6. 
The overall time resolution is given by 

At rcs = max(At rcs , At dctcctor ) . (18) 

The photon flux depends on the energy interval and therefore it is sensitive 
to the spectrum of the burst. A good phenomeno logical fit for the high energy 
GRB photon spectrum is 

N{£) = Ri £~ p , (19) 

where £ is the energy in the rest frame of the burst 1 , and N(£) has units of 
photons keV" 1 s -1 . This fit is valid for energies higher than £ , where £q is 

1 The notation in this section is the following: 8 denotes energies in the rest frame of 
the burst and E energies measured in the Earth frame. 
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typically of order of a few hundred keV. In what follows, we will take £q 
100 keV. For this spectrum the flux in the energy band (Ei, E 2 ) is 

Ri f E ^ l+z ) n _ p 



PJE 1 ,E 2 ) = — - — — / £~ (j d£. (20) 

7V 11 Aird(z) 2 1 + z Je 1 (i+ z ) V ; 

The factors (1 + z) in the limits of the integral transform rest-frame energies 
into Earth measured energies, and there is an extra (1 + z) because of the 
cosmological time dilation; d{z) is the cosmological distance. 

Let us introduce the (isotropic equivalent) peak luminosity in the GRB 
rest frame is _ 

roc 

L veak = R x \ £~' 3 £d£. (21) 



I So 

Combining Eq. l2T?land Eq. |^we obtain 

Atres _ 5 4vr^) 2 /3-l {l + zf _ ^ 

A ft 2 C-peak E\ ^ — E\ 

In the following sections we will also need to estimate the resolution from 
the photon energy flux P(Ei,E 2 ) and from the luminosity of the bursts. 
Using 

P(E 1 ,E 2 ) = 1 / £-f£d£, (23) 

4ird(z) 2 (1 + z) Je 1 (i+z) 

the time resolution is obtained by combining Eq. I2H 1221 an d 12*3*1 

b 13-1 n , E%~ p - E\~P 
Atrcs ~ A P{E U E 2 )(3-2 [1 + Z) E\-V - E\- p ' ( } 

3.2 A bound on £ 

Comparing Eq. ED and 1221 we find that we can test symmetry breaking scales 
up to: 

AH (3 - 2 p _ 2 f fE 2 \^\ Afi n 
* < 8^ J=l 8 ° Cpeak i 1 " fej ) eT Gn(Z) ' (25) 

where all the redshift dependence is contained in the function Q n {z) (remem- 
ber that Cpeak is defined in the rest frame of the burst). 
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Figure 2: Q n {z) as a function of the redshift for (3 = 2.5. For small redshifts 
Q n {z) has two branches, the upper one corresponds to Ai detector = 10~ 4 s and 
the lower one corresponds to At detector = 10~ 3 s. 



For distant bursts, the limiting resolution is determined by the photon 
arrival rate and therefore At res = At res . In this case, the function Q n {z) is 
given by: 



Gn(z) 



[l + n) 



n m (i + zf + n A 



[l + z) n dz 



Vfim + ^Jn m (i + zf + n A 

dz 

(1 + z) 



(26) 



On the other hand, for nearby bursts At res is limited by the detector res- 
olution. Fig. |21 depicts the behavior of Q n {z). For small redshifts, Q n {z) 
increases from zero up to a maximum value where Ai detector = At res . Q n {z) is 
shown in this regime for Ai de tector — 10~ 4 s (upper branch) and 10~ 3 s (lower 
branch). For higher redshifts, At res is dominated by At res . The function 
Q n {z) decreases then up to a minimum and finally increases again at high 
redshift (for n — 1 the growth begins at z > 5, which is too large to be 
experimentally interesting). As the limit on £ is proportional to Q\l n {z) the 
best limit is obtained at small redshifts, when At res = At de tector- 

This is a non intuitive result. One would expect that the best bound on 
£ is obtained from photons arriving from bursts at very high redshifts, which 
have significant larger delays than photons from nearby bursts. However the 
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distance also dilutes the photons decreasing the time resolution. Moreover, 
due to the redshift of the energy, a fixed energy band on Earth corresponds to 
an intrinsically higher energy band which is more scarce in photons. These 
two combined effects overcome the improvement of the bound due to the 
larger delay. Thus, to obtain the tightest bound it is preferable to use low 
redshift bursts. As it can be seen in Fig. |2l low redshift means here z < 1 
(but this depends, of course, on the intrinsic detector resolution). Note that 
we have not yet taken into account the possible attenuation of the interstellar 
medium. This result is merely based on geometric considerations. 

Let us turn our attention to the influence of the energy range in which we 
observe the burst. If we observe photons in the interval (Ei, E 2 ), fi n satisfies 



The exponent n + 1 — (3 is positive for n > (3 — 1. The parameter (3 can vary 
between 1.6 < f3 < 5 [221 (^ ne lower limit of (3 ~ 5 corresponds, however, 
to bursts which do not have a high energy tail and are consequently not 
interesting for our purposes). For typical bursts (3 ~ 2.5. For n > 2 and 
(3 ~ 2.5 the exponent n + 1 — (3 is positive and therefore the best bound is 
obtained by observing in the highest possible energy band. On the contrary, 
for n — 1 and (3 ~ 2.5 the exponent is negative and it is advantageous to 
use the low energy bands. Note that this analysis is based on the assumption 
that the spectrum is given by Eq. EH which is valid for energies higher than 
Sq. We should therefore always observe at energies higher than £ Q . 

Like the discussion about the optimal redshift, this is another non intu- 
itive result. Contrary to what would be naively expected, we have shown that 
the optimal energy range is not necessarily the highest one, but depends on 
the model of symmetry breaking and on the burst spectrum. These results 
apply only for a given detector with a fixed collecting area. Observations 
in a higher energy band might be advantageous under different conditions, 
for instance, if they are made with a different telescope with a larger area. 
Finally, cosmic attenuation sets an upper limit on the energy range. We 
discuss this issue in the section EHfl 




(27) 



and the energy dependence in Eq. is bounded by 




(28) 
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From Eq. [2E1 is evident that the larger the ratio E 2 /Ei, the better the 
bound (broadening the energy band increases the number of photons and 
consequently the telescope resolution); however, in order to avoid a large 
spread in the arrival times, the detected photons in each channel should 
have comparable energies. 

The bound can be rewritten as 



l/n 



£>peak n i \ 
Qn{z) 



where 

AH {3-2 p_ 
>&m£c2 13-1 * 




(29) 



l/n 

(30) 



We have introduced the quantity £* = 6.3 ■ 10 51 erg s _1 , which will be useful 
later on when dealing with luminosity distributions. In Eq. 1221 we have 
explicitly separated the redshift and burst luminosity and included all the 
numerical values and telescope dependent quantities in the constant a, 

fa A x l/n / 771 N / , 771 N l-/3\ 1 / n 



a ~ 10 U \(3-l) UoOOcmV VMevJ ^ \E 1 

(31) 

We see again that for n = 1 and (3 > 2 going to high energies does not 
improve the bound. 

The order of magnitude of achievable bounds for typical bursts obtained 
in Eq. H2 is 10~ 3 for n = 1 and 3 • 1(T 13 for n = 2 (for (3 = 2.5 in both 
cases). This is in agreement with actual limits found for specific bursts in 
the literature. Ellis et al. ^2] used a wavelet analysis to look for correlations 
between redshift and spectral time lags between the arrival times of flares 
at different energies, and obtained the bounds of £i > 5.6 ■ 10 -4 and £2 > 
2.4- 10 -13 at a 95% of confiance level. Subsequently this result was improved 
[Tlj by using a more complete data set of transient sources with a broad 
spread in redshifts to correct for intrinsic time delays (we discuss intrinsic 
delays in section |3j). A different approach was adopted by Boggs et al. [T7| . 
who considered a single and extremely bright burst, GRB 011206. This burst 
yielded the bounds of £1 > 0.015 and £ 2 > 4.5 • 10~ 12 respectively. 
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Figure 3: Optical thickness of the intergalactic medium at 10 GeV 

3.3 Cosmic attenuation 

At energies of TeV and higher, the universe becomes opaque due to the 
interaction of the gamma ray photons with the background light to create 
electron-positron pairs, 77 — > e~e + . The cross section of this reaction is 
maximized when the product of the energies of both photons is ~ (m e c 2 ) 2 . A 
photon of 10 TeV will interact with an infrared photon, for instance, creating 
an electron-positron pair. 

For n = 1 we are concerned with lower energies, from a few hundreds 
keV to a few MeV. A these energies we can safely neglect attenuation. Fig El 
shows the optical thickness of the extragalactic medium at 10 GeV [20J. As 
it can be seen, attenuation only becomes important at high redshift, z > 3. 
Clearly, it is safe to ignore this effect at lower energies. 

3.4 A quantitative example : GRB 050603A 

For an observed burst on Earth, we can skip Eq. ESI which depends on 
quantities defined in the rest frame of the burst (C pea k and £q) and calculate 
the bound directly from the flux on Earth. The quantities usually measured 
are the energy flux P(E±, E2) (in erg cm -2 s" 1 ) or the photon flux P 1 (E\, E2) 
(in photons cm -2 s _1 ). 

GRB 050603A is a very bright burst with measured redshift z = 2.821 
observed with two satellites : Konus-Wind (area ~ 200 cm 2 ) and Swift-B AT 
(area ~ 5200 cm 2 ). The time-integrated spectrum is well fitted by a high 
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energy photon index (3 = 2.15. Konus measured a peak flux of 3.2 • 10 5 erg 
cm -2 s _1 in the (20 keV, 3 MeV) band. These values correspond to: 



At rcs (20 keV, 3 MeV) 


~ 5 ■ 10" 4 s 


n — 1 


At del (20 keV) 


~ 9 ■ 10- 6 /£ s 


n = 2 


At del (20 keV) 


~ 5 • io- 29 /e 2 s 



For the same burst Swift detected a peak flux of 31.8 photons cm 2 s 1 in 
the (15 keV, 350 keV) band. This implies 



At rcs (15 keV, 350 keV) 


~ 4 ■ 10" 5 s 


n — 1 


At del (15 keV) 


~ 7-10- 6 /^ s 


n = 2 


At del (15 keV) 


~ 3 • io- 29 /e 2 s 



If the time delays at these low energy bands were not contaminated with in- 
trinsic spectral delays at the source, this burst could have led to the following 
limits: 





Konus 


Swift 


n = 1 


f > 0.02 


£>0.2 


n = 2 


£ > 3- 10~ 13 


£ > 9 ■ 10~ 13 



The time resolution used in this example, 4 ■ 10 -5 s, is in fact below the 
limiting time resolution of Swift, 0.1 ms. Swift resolution would have led to 
a slightly weaker bound, £i > 0.07. Furthermore, we have considered here 
photons all the way down to 15 keV (Swift) or 20 keV (Konus). Limiting the 
discussion to photons above ~ 100 keV, for which Eq. [THlis accurate, would 
have reduced the limits further. 

The numbers obtained here are clearly idealized and should serve only as 
an example of what can be achieved. In order to set an experimental bound, 
we need to observe the burst in two different energy channels and compare 
the arrival time in each one. An additional problem is posed by the observed 
intrinsic delay in the emission times of the photons at the source 121] 
which we discuss in section El 
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Figure 4: The luminosity and the redshift needed to give a bound of £o = 0.2 
for n — 1. Higher bounds are obtained above the curve. 



4 Distribution of bursts 

GRB 050603A could have given a very powerful bound on £. But how likely 
would it be to detect a burst yielding such a bound or a higher one ? From 
the results of the previous sections, the more luminous and closer the burst, 
the stronger the bound. In this section we estimate the probability of finding 
such a burst, given an empirical luminosity and space distribution of bursts. 

To improve a bound, £o, we need to detect a burst with a luminosity and 
a redshift such that 

° n £ Qn{z) > Q . (32) 

Eq. E21 define a region in the luminosity and redshift space-phase (see Fig. 
EJ). The probability for such a burst to happen is given by the integral of the 
probability of finding a burst over this region. 

Let us introduce the local peak luminosity function, $o(F), defined as the 
fraction of GRBs with luminosities in the interval logL and logL + dlogL, 
can be approximated by 



$o(L)-c j {L/ur u<h<AiU , (33) 
where Cq is a normalization constant such that the integral over the luminosity 
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function equals the unity. 

There is a strong evidence that long GRBs follow the comoving star for- 
mation rate (SFR), Rsfr(z). Namely, Rgrb(z), the comoving GRB rate 
satisfies Rgrb(z) oc Rsfr(z). We employ the parametrization of Porciani et 
al. JH] for the comoving SFR distribution. From it we write 



23exp(3.4z) yUn(l + zf + ^ A 
RgRb{z) ~ P0 exp(3.4^) + 22 (TT7F ' ( 3 

The luminosity function at redshift z is therefore $(z, L) = $o(L) Rgrb(z). 
Guetta et al. j2S] used the BATSE peak flux distribution, to estimate the 
parameters po, «i and «2- For long bursts they found two different fits: 





«2 


Ax 


A 2 


(erg s _1 ) 


Po (Gpc~ 3 yr" 1 ) 


I -0.1 


-2 


30 


50 


6.3 • 10 51 


0.18 


II -0.6 


-3 


30 


50 


1.6 • 10 52 


0.16 



Short GRBs, which constitute about one quarter of the observed bursts, do 
not follow the SFR fZE\ 12*7] and will not be discussed here. 

The probability of detecting a burst which sets a bound £ > £o is 

m > « = f ^M^M* r * (L) d log L, (35) 

JO 1 + 2 dz J A 

where the factor (1 + z)~ l accounts for the cosmological time dilation and 
V(z) is the comoving volume. The factor A is defined as 

A = max ( ^°J^* . , Crain) , (36) 

where £ m i n is the minimal luminosity for a burst to be detected 

2-/3 

£ mm = Air d{zf (1 + zf til -J^L M min . (37) 

C m in depends on the sensitivity of the telescope, J\f m in, which is the minimal 
photon flux necessary to trigger the instrument. We considered for n = 1 an 
idealized detector, Det. I, with an area of 5200 cm 2 making observations in 
the energy band (500 keV, 2 MeV). For n = 2, we focused on the forthcoming 
spatial observatory GLAST which will be launched in 2006. We considered 
the energy band (100 MeV, 1 GeV) where GLAST is expected to have an 
effective area of 8000 cm 2 . 

Integrating numerically Eq. |HSl w e find 
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GRB 021206 



Det. I 


GLAST 


6 


rate (bursts yr x ) 


6 


rate (bursts yr x ) 


0.015 


(7.4 - Q.9)-10-' 


4.5 ■ 10 12 


9.4 - 12.8 


0.05 


(2.8 - 3.1)-10 1 


7 ■ lO" 12 


1.9 - 2.2 


0.1 


(9.1 - 6.9)-10^ 2 


io- n 


(4.8 - 3.5)-10- x 


0.5 


(2.1 - 1.3)-10- 3 


2.5 • lO" 11 


(5.0 - 2.4)-10~ 3 


1.0 


(3.5 - 2.0)-10" 4 


5- lO^ 11 


(7.2 - 3.9)-10- 5 



The values in the first row correspond to the bounds obtained from GRB 
021206 jTTJ . To estimate C min) we assumed that the ideal detector, Det. 
I, has a sensitivity of M m i n ~ 1 ph cm -2 s" 1 , which is comparable to the 



sensitivities estimated by Band [28] for several detectors in similar energy 
bands. For GLAST, we took M m in ~ 4 ■ 10 -6 ph cm -2 s _1 , which roughly 
means that the detector is very quiet at these energies and a detection of 6 
photons during the whole duration of the burst (~ 3 minutes) is enough to 
identify it. 

Finally we took into account the partial sky coverage of any real telescope 
to compute the number of bursts observed per year. GLAST opening angle 
will be ~ 2 stereoradian; we used a similar opening angle for our idealized 
detector. 



5 Comparison with other works 

The idea of using GRBs to set experimental bounds on a possible violation of 
Lorentz symmetry was first suggested by Amelino-Camelia et al. jH]- Later 
on, several groups made use of it to explore these limits [T31 El Ej. The 
current best bounds have been obtained by Boggs et al. [T7j who used a very 
bright burst, GRB 021206, to set a limit on the symmetry breaking scale. 
The data consisted of light curves in six energy bands spanning 0.2 - 17 
MeV. The redshift of this burst is not known, but an approximated redshift 
of z ~ 0.3 was estimated from the spectral and temporal properties of the 
burst (this method involves, however, a very high uncertainty which can be 
as high factor of 2). 

The observed fluence of GRB 021206 is 1.6 ■ 10 -4 ergs cm 2 at the energy 
range of 25-100 kev [21]. This puts GRB 021206 as one of the most powerful 
bursts ever observed. GRB 021206 also shows a very atypical photon spec- 
trum at the MeV range. Instead of decreasing with the energy following a 
power law with (3 ~ 2.5, it is almost flat from 1 MeV up to 17 MeV, namely 
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(3 ~ (this implies in particular that F v increases with v in these energies). 
This flatness allows to resolve a fast flare and to determine its peak time 
and uncertainty in several bands. The analysis of the dispersion of these 
peak times yields to the lower bounds £1 > 0.015 and £2 > 5 ■ 10~ 12 [TT] . 
Applying our method on the energy band 15-350 keV using data from the 
GRB 050603A, we obtained a theoretical upper limit to the lower bounds 
of £1 > 0.2 and £2 > 9 ■ 10 -13 . These numbers represent the best bounds 
that could be obtained if the time resolution of the detector was high enough 
(~ 5 ■ 10~ 5 s), the detector noise was negligible and we had at our disposal 
the light curves in at least two energy channels. 

Our conclusions on the optimal redshift and energy band are based on a 
power law spectrum E~~P with (3 > 2. They arise from comparing the time 
delay, which always increases with the energy, with the time resolution of 
the telescope. Comparing both energy dependencies, we found in section I3~2l 
that for n — 1 and (3 > 2 is better to observe at low energies. However if 
(3 < 2, like in GRB 021206 in the MeV range, this conclusion does not hold 
and it is preferable to use the highest available energy band. In this case the 
GLAST Large Area Telescope will be a very powerful tool. It is expected 
to be sensitive from 20 MeV to 300 GeV with a peak effective area in the 
range 1-10 GeV of 8000 cm 2 . Observing with GLAST in energy bands below 
10 GeV where cosmic extinction is still negligible (see fig. EJ) can improve 
dramatically our current bounds. At present little is known about GRB 
emission at energies higher than 50 MeV and therefore it is not possible to 
estimate how common are bursts with (3 <2. Taking advantage of atypical 
bursts to explore even higher energies is an exciting possibility to keep in 
mind, however at present it is difficult to design a strategy based only in 
these bursts. 

As already mentioned, our bound must be interpreted as a theoretical one, 
i.e. the highest bound that could be set, were the best conditions achieved. 
We already commented on the necessity of observing in at least two energy 
channels and to take into account the real time resolution of the detector. 

An additional serious problem is the intrinsic lack of simultaneity in the 
pulse emission in the keV regime ^3 121] • Soft emission has a time delay 
relative to high energy emission [22] • While the reason for this phenomenon is 
not understood, an anti-correlation between the spectral evolution timescale 
and the peak luminosity has been found [T5] . There are two different ways 
to deal with the intrinsic delay. The first is to try to reduce it by choosing 
very luminous bursts and observing in MeV or higher, where the delay, if still 
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exists, seems to be smaller. The second approach is based on the fact that the 
delays produced by a violation of Lorentz symmetry increase with the redshift 
of the source, whereas intrinsic time delays are independent of the redshift 
of the source [T3j. Thus, a systematic comparison of a delays in a group of 
bursts with known redshifts could enable us to distinguish between intrinsic 
and redshift dependent delays. Using a sample of 35 bursts with known 
redshifts, Ellis et al. [Hj established a lower limit of £1 > 7 • 10 -4 on the 
symmetry breaking scale. These bounds are two orders of magnitude lower 
than our theoretical limits. This difference demonstrates the importance that 
intrinsic time delays, noise and the real instrumental resolution can have. 

6 Conclusions 

Our goal was to explore the potential of GRBs to set bounds on Lorentz 
violation and to find optimal techniques to do so. We modified the dispersion 
relation of photons by adding an extra term proportional to the photon 
momentum to the power n + 2. We have shown that in models with n — 1 
it is possible to explore energies which are close to the Planck energy. When 
n — 2, the energies explored are smaller, around 10 7 GeV. These bounds 
are idealized and they do not take into account experimental limitations 
or the intrinsic time-structure of the 7-ray emission. They should serve as 
theoretical estimations of what can be achieved. The methodology we use 
here can be used to design future optimal experiments for observing this 
effect (or setting bounds on it). 

We have modelled the burst high energy emission with a power law spec- 
trum E~P with (3 > 2. This fits well the time integrated emission of most of 
the bursts. We found two non intuitive results: (i) The optimal redshift to 
set the strongest bound is less than 1. (ii) For n — 1, low energy, rather than 
high energy emission is preferred. Both results are counter-intuitive since 
the Lorentz violation delay increases with the distance and with the energy. 
However, distance or observations at high energies (where the flux is lower) 
dilute the photons reducing the temporal resolution achieved on Earth. It 
turns out that this is the dominant effect. 

In the models with n — 2, going to higher energies always improves the 
bounds. Here the situation will be remarkably changed when the spatial 
observatory GLAST will become operational. 

We have also investigated the probability of improving the current exper- 
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imental bounds, given a phenomenological luminosity and space distribution 
of bursts. As we are discussing idealized bounds, this probability should only 
be trusted up to an order of magnitude. 
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